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Nonthermal Emissions from Particles Accelerated by Turbulence 

in Clusters of Galaxies 
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ABSTRACT 

We consider nonthermal emission from clusters of galaxies produced by parti- 
cle acceleration by resonant scattering of Alfven waves driven by fluid turbulence 
through the Lighthill mechanism in the intracluster medium. We assume that 
the turbulence is driven by cluster mergers. We find that the resonant Alfven 
waves can accelerate electrons up to 7 ~ 10^ through resonant scattering. We 
also find that the turbulent resonant acceleration can give enough energy to elec- 
trons to produce the observed diffuse radio relic emission from clusters if the 
clusters have a pool of electrons with 7 ~ 10'^. This mechanism can also ex- 
plain the observed hard X-ray emission from clusters if the magnetic field in a 
cluster is small enough (< fiG) or the fluid turbulence spectrum is flatter than 
the Kolmogorov law. The fluid turbulence could be observed with Astro-E2 in 
the regions where diffuse radio emission is observed. Although non-gravitational 
heating before cluster formation (preheating) steepens a relation between radio 
luminosity and X-ray temperature, our predicted relation is still flatter than the 
observed one. 
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1. Introduction 

The intracluster medium (ICM), the largest detected baryonic component of clusters of 
galaxies, contains thermal gas characterized by temperatures in the range of ~ 2 — 10 keV 
and by central densities of ~ 10^^ cm^'^. In addition to the thermal plasma, observations 
have showed that the ICM contains nonthermal elements. Diffuse synchrotron emission from 
the ICM has been observed in many clusters (e.g., Kim et al. 1990; Giovannini et al. 1993; 
Giovannini & Feretti 2000; Kempner & Sarazin 2001). Diffuse radio sources in clusters 
are often classified either as peripheral cluster radio relic sources, or central cluster radio 
halo sources. These radio halo and relic sources are not associated directly with individual 
galaxies, but appear to be produced by relativistic electrons in the intracluster space. The 
typical energy of the electrons that produce the radio emission depends on the observed 
frequency and on the intracluster magnetic field, but it is estimated to be several GeV. 
Moreover, hard X-ray emission has been detected as a nonthermal tail at energies >20 
keV in at least two clusters. The Coma cluster was detected with both BeppoSAX (Fusco- 
Femiano et al. 1999) and the Rossi X-Ray Timing Explorer (RXTE ; Rephaeli, Gruber, & 
Blanco 1999). Possible weaker excesses may have been seen in Abell 2199 (Kaastra et al. 
1999), Abell 2319 (Gruber & Rephaeh 2002), and Abell 3667 (Fusco-Femiano et al. 2001). 
This radiation is often believed to be inverse Compton emission, in this case produced by 
relativistic electrons with the energy of ~ 5 GeV and Lorentz factors 7 ~ 10"^. 

There are several possible origins for the relativistic electrons (Atoyan & Volk 2000; 
EnBlin & Gopal-Krishna 2001; EnBlin & Briiggen 2002). Electron acceleration at shocks 
in the ICM has been the most popular idea (e.g., Colafrancesco & Blasi 1998; Bykov et 
al. 2000; Totani & Kitayama 2000; Waxman & Loeb 2000; Ostrowski & Siemieniec-Oziebo 
2002). These shocks may be attributed to the interaction between jets originated from active 
galactic nuclei (AGNs) and ICM, to galactic winds, or to the accretion of intergalactic gas. 
Since most of the observed radio sources are very large and are often found in irregular 
clusters, they may be attributed to shocks formed by cluster mergers (Roettiger, Stone, 
& Burns 1999; Blasi 2001; Buote 2001; Fujita & Sarazin 2001). Takizawa & Naito (2000) 
calculated the nonthermal emission from relativistic electrons accelerated around the shocks 
produced during a merger of clusters with equal mass. They found that the hard X-ray 
and radio emission is luminous only while signatures of merging events are clearly seen in 
the ICM. For some clusters, observations suggest that electrons are actually accelerated 
by shocks. For example, Markevitch & Vikhlinin (2001) showed that a high-temperature 
shocked region is associated with diffuse radio emission in Abell 665 by analyzing data 
from the Chandra Observatory. Shock acceleration or re-acceleration seems one of the most 
promising candidates of the source of relativistic particles in the localized, peripheral radio 
relic sources. 
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On the other hand, some clusters have large diffuse radio emission that does not seem 
to be associated with shocks. Since the lifetimes of the particles which produce the radio 
emission in clusters are short, it is unlikely that radio emission will be found far from shocks 
if particles are only accelerated at shocks. Since the large radio sources arc generally found in 
merging clusters, the simplest explanation would be an acceleration mechanism associated 
with mergers which is not directly related to shocks or which persists in the plasma for 
some period after the passage of a merger shock. For Abell 3667, a diffuse radio emission is 
detected in the northwest of the cluster center (Rottgering et al. 1997). Although this cluster 
has a rapidly moving substructure, the Mach number is estimated to be only ~ 1 (Vikhlinin, 
Markevitch, & Murray 2001). Thus, even if there are shocks in the cluster, the efficiency of 
shock acceleration may be too small to produce the observed synchrotron spectrum within 
the standard theory of shock acceleration. Moreover, although Abell 2256, which is known 
as a merging cluster, also has a large diffuse radio emission (Rottgering et al. 1994), the 
X-ray temperature map does not show shock features (Sun et al. 2002). The viewing angle 
may prevent us from observing the shock. However, even if particles arc accelerated at a 
bow shock, the size of the radio emission (~Mpc) seems to be too large for a cluster to cross 
within the lifetimes of the particles responsible for the emission. Recently, Gabici & Blasi 
(2002) investigated the effect of the low Mach numbers of merger-related shocks on particle 
acceleration at the shocks. They suggest that major mergers, which often invoked to be sites 
for the production of extended radio emissions, have shocks which are too weak to result in 
appreciable nonthermal activity. 

Fluid turbulence in ICM is another possible origin of particle acceleration (Blasi 2000; 

Brunctti et al. 2001a,b). Numerical simulations show that cluster mergers can generate 
strong turbulence in the ICM (Roettigcr et al. 1999; Ricker & Sarazin 2001). The fluid 
turbulence can induce MHD waves (Lighthill 1952; Kulsrud 1955; Kato 1968; Eilek & Hen- 
rikscn 1984), which accelerate particles through wave-particle resonance (Jaffc 1977; Roland 
1981; Schlickeiser, Sievers, & Thiemann 1987). Contrary to shock acceleration, the turbu- 
lent resonant acceleration does not require strong shocks. Even if a subcluster is moving 
with a subsonic velocity in a host cluster, the induced turbulence may produce high energy 
particles and nonthermal emission. As the first step of this study, in this paper we consider 
particle acceleration by Alfven waves, which are easily generated by fluid turbulence and 
do not damp easily. The effects of other MHD waves may also be important (Schlickeiser 
& Miller 1998; Ragot & Schlickeiser 1998). Recently, Ohno, Takizawa, & Shibata (2002) 
studied the resonant acceleration in the Coma cluster; we will study the turbulence and ac- 
celeration in clusters more generally. As will be shown in §3, turbulent resonant acceleration 
is not effective when two clusters coUide with a large relative velocity because the lifetime 
of the induced fluid turbulence is shorter than the acceleration time of the particles. Thus, 
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turbulent acceleration is most effective when two clusters approach from a relatively small 
initial separation because the relative velocity when they pass each other is then relatively 
small. Moreover, this mechanism is more likely to apply to peripheral relics rather than 
central halos because two clusters have the largest relative velocity at the point of closest 
approach of the cluster centers; this large velocity decreases the time-scale of the fluid turbu- 
lence and also makes ram-pressure stripping of the gas of the smaller cluster effective, which 
reduces the volume of the turbulent region. We will concentrate on models and parameters 
appropriate for the cases where this mechanism is likely to be most important. 

This paper is organized as follows. In §2, we summarize our models. In §3 we give the 
results of our calculations and compare them with observations. Conclusions are given in 
§4. 



2. Models 

2.1. Fluid Turbulence, Lighthill Radiation, and Particle Acceleration 

We use the turbulent resonant acceleration model studied by Eilek & Henriksen (1984); 
we summarize it in this subsection. We assume that fluid turbulence is induced by the motion 
of a smaller cluster in a larger cluster and its energy spectrum is described by a power law, 

Wf{K) = W^K-"" , (1) 

where k = 2n/l is the wavenumber corresponding to the scale /, Wf{K,)dK, is the energy per 
unit volume in turbulence with wavenumbcrs between k, and k, + dn, and Wj and m are the 
constants. If one expresses the turbulent spectrum in terms of eddy size, the spectrum is 
represented by Wf{l) oc Z'""^. The cascade of the fluid turbulence extends from a largest 
eddy size Iq — 2ti / kq down to a smallest scale determined by dissipation. Id ~ lo{Re)~^f^ , 
where Re is the Reynolds number. Since most of the energy of fluid turbulence resides in 
the largest scale, the total energy density of fluid turbulence is given by Et ~ pv^, where p 
is the fluid density and Vt is the turbulent velocity of the largest scale Iq. The normalization 
can be derived from the relation Ef — J Wf{K)dK and is 

W^^EtK^p-'/R, (2) 

where 

R= ^ (3) 

m — 1 KTWf[KT) 

In the above equations, kt — 2t:/ It and It is the wavelength below which Alfven waves are 
driven. 
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Fluid turbulence will generate Alfven waves via Lighthill radiation (Lighthill 1952; Kul- 
srud 1955; Kato 1968). The maximum wavelength of Alfven waves driven by the turbulence 
cannot be larger than Iq. Instead, Alfven waves are expected to be driven at the wavelength 
at which there is a transition from large-scale ordered turbulence to small-scale disordered 
motions. Unfortunately, there is no absolute definition of this transition. Following Eilek & 
Henriksen (1984), we adopt the Taylor length as an estimate of the transition scale. The 
Taylor length is given by 



ll 



dx 



lo{lh/Re 



,1/2 



(4) 



where Vf is the turbulent velocity vector, x is the space coordinate, and Vf^^ is the x- 
component of the vector Vf (Tennekes & Lumley 1972; Eilek k, Henriksen 1984). The Taylor 
scale is determined by the ratio of the r.m.s. turbulent velocity {[vj^^^'^)^ which is determined 

by the large scale motions, to the r.m.s. velocity shear {[{dvf^x/d^TY^'^)-, which is mainly 
due to small scale motions. 



A fluid eddy of size I has a velocity Vf 

Vf{l) 





1/2 


. p . 





|Vf|, 

pR 



K 

Kt 



l—m 



1/2 



Turbulence on a scale k will radiate Alfven waves at the wavenumber 



(5) 



(6) 



where va = B/^/A7^p is the Alfven velocity in plasma with the magnetic field B and the gas 
density p. Let Vf[K{k)] be the fiuid velocity on the fiuid scale K{k) which drives Alfven waves 
of wavenumber k. From equations (5) and (6), we obtain 

/ \ l/{3-m) / ^ \ (l-m)/(3-m) 

\pv\Rj \kt 



Vf[K{k)] ^VA 



(7) 



We assume that the energy going into Alfven waves at wavenumber k is given by a 
power law 

lA(k) = hik/KT^ , (8) 

where lA{k)dk is the energy per unit volume per unit time going into Alfven waves with 
wavenumbers in the range k ^ k + dk, and Iq and St are the constants. In this case, the 
total power going into the Alfven mode from fluid turbulence is 

Pa - r^lA{k)dk^^(^Y\ (9) 

Jkr 



St 
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where kr — [vf{lT)/vA]i^T- We have assumed that kr <S k^ax and St > 1. On the other 
hand, according to the Lighthill theory, Pa is given by 

Pa - Va"^^ X fWfihfKT , (10) 
Va 

where tja is an efficiency factor of order unity (Kato 1968; Henriksen ct al. 1982; Eilek & 
Henriksen 1984). By comparing equations (9) and (10) and using equations (6) and (7), it 
can be shown that 

= . (11) 



3 — m 



h = VA{st - l)fyv\ -4- , (12) 



and 

f E \^ 

(Eilek & Henriksen 1984). 

The Alfven waves respond to this energy input according to 

dWA{k) 



dt 



lA{k)-'y{k)WA{k), (14) 



where WA{k)dk is the energy per unit volume in Alfven waves with wavenumbers in the 
range k ^ k + dk. and ^{k) is the damping rate (Eilek 1979; Eilek & Henriksen 1984). We 
consider damping by the acceleration of relativistic particles. For a particle with momentum 
p, the resonance condition is given by 

k ^ Q.me/[p{iJ,- va/c)] , (15) 

where Q = eB/rrieC is the nonrelativistic cyclotron frequency, and fi = v\\/c is the projected 
particle velocity along the magnetic field normalized by the light velocity c (Eilek 1979). 
Then, the damping rate due to accelerating particles is 



7acc(A:) = / p2 



1 _ ( !ji + ^rue ^ ^ 
; pk 



-^dp, (16) 



where e and rrie are the electron charge and mass, and Prain{k) = flme/[k{l — va/c)] from the 
resonance condition (eq. [15]). The electron phase space distribution function is f{p) (Eilek 
1979; Eilek & Henriksen 1984). In general, pmax ^ Pmin (see equation [40] below) and the 
electron phase space distribution function declines rapidly at large p, so the upper limit does 
not contribute significantly to the integral in equation (16). In fact, although the integral 
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converges whenever /(p) is steeper than p ^, we later consider only the case of 4 < s < 5, 
where f{p) oc p~^. 

We consider the case where the wave and particle spectra are represented by power-law 
functions, namely, 

WA{k) = W^k-" , (17) 
f{p) = hp-' , (18) 

where W\ and /o arc the constants. Since the characteristic wave response time is ~ 7^;^^, 
which is short (~ 10^ yr) compared to other relevant time-scales in our model, we can assume 
that dWA/dt — 0. Thus, from equations (14) and (16), and assuming va <S c, we obtain 



fo ^ 87r3 e2^;2 
1/ — s + St — 3 . 



s-2 



(19) 
(20) 



On the other hand, the electron distribution evolves according to the equation. 



df _ 1 d 
dt dp 



p'Dip)^ + Sp'f 



(21) 



Here, S = 4(5^ + i?c]y[g)e^/(9mgC^) is the inverse Compton and synchrotron loss coefficient, 
and -BcMB = 3.25(1 + fiG. In this paper, we assume that the redshift of a model cluster 
is 2; = 0. The diffusion coefficient owing to the resonant scattering off random Alfven waves 
is 



D{p)^ 



WA{k) 



(p) 



k 



Va flm^ 
c pk 



dk , 



(22) 



where A;min(p) = ^me/lpil — va/c)]. The value of the upper cutoff /Cmax is not important as 
long as kmax ^ ^min and u > 0. Using equations (17) and (18), equation (21) is rewritten as 



9/ _ ]_d_ 
dt p^ dp 



Ap'^^% + Sp'f 
op 



with 



A 



S-2-V 



2t:v{v + 2) c 
Equation (23) can further be rewritten as 

df _ f{p) f{p) 

dt ta{p) te{p) 



fo 



T ■ 



(23) 



(24) 



(25) 
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with 

ta{p)^[s{s-u-l)Ap^-Y\ (26) 
Up)^[(s-A)Sp]-' . (27) 



If i/ = 3, equations (26) and (27) show that the ratio 

tgjp) _ 1 S 
te{p) S A 



(28) 



is independent of p. In other words, if and only ii u — 3, the form of /(p) does not evolve 
(Eilek & Henriksen 1984). Moreover, Eilek & Henriksen (1984) showed that even ii u ^ 3 
initially, Wa and / evolve so that u and 5"/ (sA) approach three and unity, respectively, on a 
timescale of approximately tp = min[ta, te]. Note that Eilek & Henriksen (1984) showed that 
equation (23) has another time-indcpcndcnt solution with f(p) oc p~^, in addition to the 
time-independent sohition with f{p) oc p^*' (s = S/A) which we adopt here. However, the 
former solution requires externally supplied fluxes of particles in contrast with the latter. 
Since we do not consider external sources of particles in this paper, we ignore this other 
solution. 

Prom equations (6) and (7), we obtain 
Since the sum of synchrotron and inverse Compton emission rate is given by 

/"Pmax 

Psi= / ^7rSp^cf{p)dp, (30) 



it can be shown that 



Psi 2s(s-2) f kxP^i^A^ ' S 



Pa 15 V ^rrip J sA 



(31) 



from equations (13) and (29) when i> = 3 (Eilek & Henriksen 1984). From the resonance 
condition (eq. [15]), 

Pmax ~ flme/kr ■ (32) 
Moreover, as mentioned above, S/ (sA) a; 1 in a steady-state. Thus, in this case, 

Psi _ 2s{s - 2) 



Pa 15 



(33) 



(Eilek & Henriksen 1984). We note that equations (10) and (33) indicate that the emission 
power, Psi, mainly depends on the energy injection at the scale of It] it does not much 
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depend on the particle spectral index s, which is expected to be of the order of 4. This 
means that even if the fluid spectrum index m for k > kt is different from that for k < kt 
because of the back reaction of particle acceleration (Ohno et al. 2002), the luminosity of 
nonthermal emission does not change significantly, although the emission spectrum may be 
affected. 



2.2. Viscosity of the ICM 

In order to apply the turbulence acceleration model in §2.1, wc need to find the Reynolds 
number of the ICM, Re = loVt/fRi where vk is the kinetic viscosity. The viscosity in turn 
is represented by = ""p^cfr/S, where Up is the thermal velocity of protons and Acfr is the 
effective mean free path of protons. For the transverse drift of protons in magnetic fields, 
the mean free path is given by Aeg = A^/Ac (Spitzer 1962; Ruzmaikin, Sokolov, & Shukurov 
1989), where A^ is the proton gyroradius and Ac is the mean free path due to Coulomb 
collisions: 

A, = ^ , (34) 

'^^"47rV2n,e2 1nA ^"^^^ 
(Sarazin 1986). In the above equations, rUp is the proton mass, Ub is the Boltzmann constant, 
Tgas is the gas temperature, Up is the proton density, and In A is the Coulomb logarithm. 
Thus, we have 

= 1.3 X 10= cm^ ( I'' \ ('^) " . (36) 

Vio-3cm-3y vio^Ky \aq)\i^g) ^ ' 

We use this viscosity in §3. 

We would like to point out the relation between the viscosity and the electron maximum 
energy. As equation (32) shows, the maximum energy of accelerated electrons, Pmax, is related 
to It- Thus, Pmax should constrained by the viscosity vk or the Reynolds number Re (eq. 
Prom equations (5), (6), and (32), one can show that 

27rmec2/o"^™"^^^^t't7n ^ 



h 



/max 



BevA 



(37) 



where the maximum Lorentz factor is 7niax = Pmax/ (jTigc) for Pmax ^ nieC. Using equation (4) 
and Re — Vtlo/vK, we obtain 

- eE^Zo ^ 15^^ V^-)/^ 
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where the last two equations assume Kolmogorov turbulence (m = 5/3). If nonthermal 
radio and hard X-ray emission from clusters is due to turbulent acceleration, this upper- 
limit may explain the discrepancy between the cluster magnetic fields derived from Faraday 
rotation (e.g., Lawler & Dennison 1982; Kim et al. 1990; Goldshmidt & Rephaeli 1993; Clarke, 
Kronberg, & Bohringer 2001) and those derived from hard X-ray emission assuming a power 
law particle spectrum without an upper limit (Fusco-Femiano et al. 1999, 2000). Previous 
work indicated that a high-energy cutoff (7max ~ 10^) in the electron energy distribution 
could resolve this discrepancy (Brunetti et al. 2001a; Petrosian 2001; Fujita & Sarazin 2001). 
Our model prediction (eq. [40]) for the cutoff is close to that required, although this is a very 
rough estimate. 



2.3. Cluster Merger 

We consider a merger in which a smaller cluster falls into a larger cluster. The merger 
model is based on the galaxy infall model in Fujita (2001). 



2.3.1. Mass Profile of a Cluster 



In this subsection, we do not discriminate between the smaller and larger clusters un- 
less otherwise mentioned. The virial radius of a cluster with the virial mass My^^ and the 
formation redshift z is defined as 



3Mv 



1/3 



(41) 

_47rAc(z)pcrit(-s)J ' 

where Pcnt{z) is the critical density of the universe and Ac{z) is the ratio of the average 
density of the cluster to the critical density at redshift z. The former is given by 

Pcrit,of^o(l + 



Pcrit(^) = 



n{z) 



(42) 
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where Pcrit,o is the critical density at z — 0, and fl{z) is the cosmological density parameter. 
For the Einstein-de Sitter Universe, /^d^) — 18 tt^. 

We assume that a cluster is spherically symmetric and ignore the gravity of ICM. The 
density distribution of dark matter is represented by a power-law, 

Pm(r) = Pmv(r/rvir)"" , (43) 

where p^y and a are the constants, and r is the distance from the cluster center. The 
normalization, p^y, is given by 

3 — a , . 

Pmv = — ^AcPcrit ■ (44) 

We choose a = 2.4, because the slope is consistent with observations (Horner, Mushotzky & 
Scharf 1999). Moreover, the results of numerical simulations show that the mass distribution 
in the outer region of clusters is approximately given by equation (43) with a ~ 2.4 (Navarro, 
Prenk, & White 1996, 1997). We adopt a power-law rather than the full NFW profile to avoid 
specifying a particular value for the concentration parameter and its variation with cluster 
mass and formation redshift. At present, there is no complete and consistent understanding 
of the variation of the concentration parameter based on numerical simulations (Navarro 
et al. 1997; Bullock et al. 2001). Moreover, the central cusp that is a characteristic of the 
NFW profile does not have a significant impact on the motion of a radially infalling object 
in the cluster (Fig. 2 in Fujita 1998). For example, if we assume that the cluster mass within 
0.5 Mpc from the center is 2.5 x 10^^ Mq and that a smaller cluster falls into the center 
from an initial radius of r = 2 Mpc, the velocity difference at r = 0.1 Mpc is less than 10% 
between an NFW profile with a scale radius of 0.9 Mpc and our adopted profile with a — 2.4. 

We consider two ICM mass distributions. One follows equation (43) except for the 
normalization and the core structure; 

[1 + (r/re)2]-"/2 

PiCM(r) = p^cMMr^^^^^^^j^^^,^.^/, ■ (45) 

The ICM mass within the virial radius of a cluster is 

MicM= / 47rrViCM (r) o?r . (46) 
Jo 

The normalization picM,vir is determined by the relation = MicM/^vir, where is the gas 
or baryon fraction of the universe. This distribution corresponds to the case where the ICM 
is in pressure equilibrium with the gravity of the cluster and is not heated by anything other 
than the gravity. We call this distribution the 'non- heated ICM distribution'. We introduce 
the core structure to avoid the divergence of gas density at r = and assume Tc = O.lryir- 
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For the larger cluster into which a smaller cluster falls, we will modify the core radius to take 
account of the finite size of the smaller cluster (see §2.3.2). We use fb = 0.25(/?y0.5)~^/^, 
where the present value of the Hubble constant is written as Hq = 100 h km s~^ Mpc~^. 
The value of fb is the observed ICM mass fraction of high-temperature clusters (Mohr, 
Mathiesen, & Evrard 1999; Ettori & Fabian 1999; Arnaud & Evrard 1999), for which the 
effect of non-gravitational heating is expected to be small. 

However, X-ray observations suggest that the ICM is also heated non-gravitationally at 
least for nearby clusters (e.g., David, Jones, & Forman 1996; Ponman, Cannon, & Navarro 
1999; Lloyd-Davics, Ponman, & Cannon 2000; Xuc & Wu 2000). Thus, wc also model the 
distribution of the heated ICM using the observed parameters of nearby clusters as follows. 
In this paper, we assume that the ICM had been heated before being accreted by clusters. 
However, the distribution will qualitatively be the same even if the ICM is heated after 
accretion (see Loewenstein 2000). 

Following Balogh, Babul, & Patton (1999), we define the adiabat Kq = P/ p^"^, where 
P is the gas pressure, p is its density, and 7ad = 5/3 is the constant adiabatic index. If ICM 
had already been heated before accreted by a cluster, the entropy prevents the gas from 
collapsing into the cluster with dark matter. In this case, the ICM fraction of the cluster is 
given by 



where X34 = 10^^ erg g~^/^ cm^ (Balogh et al. 1999). 
The virial temperature of a cluster is given by 

prrip 2 Tvir 

where p — 0.61 is the mean molecular weight, and G is the gravitational constant. When the 
virial temperature of a cluster is much larger than that of the gas accreted by the cluster, 
a shock forms near the virial radius of the cluster (Takizawa & Mineshige 1998; Cavaliere, 
Menci, & Tozzi 1998). The temperature of the postshock gas (T2) is related to that of the 
preshock gas (Ti) and is approximately given by 



(Cavaliere et al. 1998). Wc assume that the gas temperature does not change very much for 
r < Tvir (see Takizawa & Mineshige 1998). In this case, the ICM temperature of the cluster 
is given by Tjcm = T2. Since we assume that the density profile of gravitational matter is 




(47) 




(49) 
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given by equation (43) with a = 2.4, the density profile of ICM is given by 

( \ - [l + (r/re)^]-3^/^ 

where j3 = (2.4/3)Tvir/TicM (sec Bahcall & Lubin 1994). Observations suggest that Ti ~ 
0.5 — 1 keV ahhough it depends on the distribution of the gravitational matter in a cluster 
(Cavaliere et al. 1998; Fujita & Takahara 2000). We choose 3Ti/2 = 0.8 keV hereafter. The 
normalization picM,vir is determined by the relation /icm = Micu/M^h- 



2.3.2. Ram-Pressure Stripping 

We first consider a radially infalling smaller cluster accreted by a larger cluster. Since 
we study the first infall of a smaller cluster before relaxation, we will ignore the infiuencc of 
the smaller cluster on the gravitational potential of the larger cluster. (Nonthermal emission 
from a relaxing cluster will be studied in §3.4.) That is, we assume that the structure of the 
larger cluster remains intact during the infall of the smaller cluster. The exception is the 
central region of the larger cluster. Since the smaller cluster has a finite size, the gravitational 
effect of the central density peak in the larger chister will be smoothed out on the scale of 
the smaller cluster's size or the tidal radius when the smaller cluster reaches the center of 
the larger cluster. The tidal radius of the smaller cluster, fti, at the distance r from the 
larger cluster center is obtained by solving the equation. 



r 



3M,(r)J ' ^^^^ 



where Ms{r) and ML{r) are the mass profiles of the smaller and the larger cluster, respec- 
tively (Binney & Tremaine 1987). The internal velocity of the smaller cluster is 



-r., ^ , , (52) 

^ti 

and the dynamical time scale is fdyn = ftiMa. As the smaller cluster approaches the larger 
cluster center, fti decreases and the infall velocity V increases. However, when the infall time- 
scale r/V, is smaller than idyn, the tidal disruption virtually stops. We define the minimum 
radius of the smaller cluster fmin as the tidal radius when r/V = idyn is satisfied first. We give 
the core radius of the larger cluster by r£,^c — inax(0.1rvir, fmin) in equations (45) and (50). 

We investigate two cases for the initial position of the smaller cluster. One is that the 
smaller cluster starts to move at the turnaround radius of the larger cluster rta- In this case. 
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we give the initial velocity of the smaller cluster, Vi, at the virial radius of the larger cluster 
r = r^yij., and it is 



2 fL,wvc fta. 

where M^^vir is the virial mass of the larger cluster. Assuming that rta = Sr^^vir on the basis 
of the virial theorem, the initial velocity is 



'GM^ . (54) 

^LjVir 

The virial radius rL,vu is given by equation (41). 

For the second case, we assume that the smaller cluster was initially located near the 
center of the larger cluster, as expected for the merger of substructure within a large cluster 
(Fujita et al. 2002). The initial distance from the center of the larger cluster is given by 
Tj = Xjrj^ vir, where Xi is a parameter and < < 1; the velocity at r = rj is zero. For both 
cases, the velocity of the smaller cluster is obtained by solving the equation of motion; 

Tt = ^ ■ ^^^^ 

for r > r^^c- We assume that Pm{T^) — Pm{fL,c) in equation (55) for r < r^^c- 

As the velocity of a smaller cluster increases, the ram-pressure from the ICM of the 
larger cluster affects the gas distribution of the smaller cluster. Sarazin (2002) indicated 
that the ram pressure affects the gas distribution at radii which satisfy 

i^ram = PL,ICM(r)i;' > Ps{f) , (56) 

where Pl,icm is the ICM density of the larger cluster, Pg is the static pressure of the smaller 
cluster, and f is the distance from the center of the smaller cluster. The static pressure 
is given by P5,icm(^^)^b^s',icm/(/^^p), where the index S indicates the smaller cluster. We 
assume that the ICM of the smaller cluster is stripped outside the radius, fgt, satisfying 
-Pram = -P^l^^st)- Howcvcr, if fti is Smaller than fgt, we reset the ICM radius to fgt = fti- 



2. 3. 3. Nonthermal Emission from Merging Clusters 

The motion of the smaller cluster relative to the ICM of the larger cluster will induce 
turbulence behind the smaller cluster. The largest eddy size and velocity are respectively 
< fst and < V, where v is the relative velocity of the smaller cluster to the larger cluster. 



- 15 - 



Thus, we define lo — fi fst and Vt — fvV, where < < 1 and < < 1. We expect 
that nonthermal emission comes from the turbulent region according to equation (33). The 
volume of the turbulent region is given approximately by 

Kurb = TTf^tVtturb , (57) 

where iturb = ^o/^t is the time-scale of the turbulence. Thus, the synchrotron luminosity and 
inverse Compton luminosity can be written as 

i'si^turb (58) 



"I" -"CMB 

and 



LJC = o2?o2 ^Sl^turb , (59) 

-°CMB 



respectively. We assume that large-scale magnetic fields are given by 

52 , 1 



^ = & X ^picuVt , (60) 



where 6 is a parameter. 



In our model, since the lifetimes of accelerated high-energy particles are equal to or 
less than the time-scale of fluid turbulence (see §3.1), the high-energy particles and the 
nonthermal emission from them are mostly confined in the turbulent region with the volume 
of Vturh behind a smaller cluster. Note that the local energy density of the high energy 
particles is smaller than the thermal energy density of ICM in the calculations in §3. 



3. Results and Discussion 

3.1. General Features 

Wc calculate the nonthermal emission from a merging cluster using the model presented 
in §2. In this subsection, we study the dependence of the nonthermal emission on the cluster 
masses. The model parameters are shown in Table 1. First, we study the case where a smaller 
cluster falls from = O.Sr^ vir, and we adapt the non- heated ICM distribution (cq. [45]) and 
the Kolmogorov law (m = 5/3) for the fluid turbulence (Model A). From now on, we assume 
Qo = 1 and /i = 0.5 for the sake of simplicity. A larger cluster forms at z = 0, and a smaller 
cluster forms at 2; = 0.5. We also assume that fi — 0.5 and b — 1. Kato (1968) showed that 
in the case of equipartition of energy and eddy size between the velocity and magnetic fields 
on small scales, rjA — 15.5/(27r). We adopt that value for tja- 
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Figure 1 shows the synchrotron luminosity of merging clusters when the smaller cluster 
is at the distance r from the center of the larger cluster. The outer bends in the curves (e.g. 
at r = 1.8 Mpc for the model of M^^^vir = 5 x IO^^Mq and M^^^ir = 1-5 x IO^^Mq) correspond 
to the radii at which the ICM radius of the smaller cluster equals to its tidal radius. When 
Miy^r = 5 X IO^^Mq and Msyn- = 1.5 x IO^^Mq, the maximum synchrotron luminosity is 
~ 3 X 10*^° erg s~^. Since the masses of the clusters are close to the maximum observed so far, 
the synchrotron luminosity should be the maximum of any clusters for the parameters we 
have assumed. This maximum luminosity is close to that obtained from observations (e.g., 
EnBlin et al. 1998; Giovannini et al. 1999; Giovannini & Feretti 2000; Kempner & Sarazin 
2001), which shows that particle acceleration originating from fluid turbulence is a promising 
candidate for the synchrotron radio emission from clusters. 

The spectral index of the nonthermal emission is given by a = (s — 3)/2. Since s = 4.5 
for m = 5/3, the spectral index for the model cluster is a = 0.75. This is somewhat smaller 
than the observed values. For instance, for 18 clusters investigated by Kempner & Sarazin 
(2001), the average spectral index is a = 1.2 at 327 MHz. One possibility is that the upper 
cutoff of the energy distribution of particles steepens the observed radio spectrum (eq. [40]). 
On the other hand, the observed steep spectrum may be affected by the emission from the 
region where fluid turbulence has almost disappeared and higher energy particles have lost 
their energy. Although it is halo emission, Giovannini et al. (1993) showed that the spectral 
index of the radio halo in the Coma cluster is flat at the halo center {a ~ 0.5) compared 
with that in the peripheral region {a ~ 1.5). In the future, more observations of the spectral 
index of halo and relic radio emission will be useful to study the turbulent acceleration there. 

Figure 1 also shows the inverse Compton luminosities. Compared to the synchrotron 
luminosities, they achieve their maxima when the smaller cluster is in the outer region of the 
larger cluster. This is because the magnetic fields in the clusters increase with the kinetic 
energy of turbulence (eq. [60]), and the fraction of the inverse Compton luminosity decreases 
as the magnetic fields increase (eqs. [58] and [59]). We present the magnetic fields produced 
by a smaller cluster at the distance r from the center of a larger cluster in Figure 2. The 
absolute inverse Compton luminosity (^ 3 x lO'^^erg s^^) we predict in Figure 1 is much 
smaller than the hard X-ray luminosities observed in several clusters. For the Coma cluster, 
which is known as a merging cluster, the BeppoSAX detected hard X-ray flux of 2.2 x 10~^^ 
ergs cm~^ s~^ in the 20 — 80 keV band. This corresponds to the luminosity of ~ 10^^ ergs 
s~^ (Fusco-Femiano et al. 1999). For Abell 2256, which is also known as a merging cluster, 
the flux of hard X-ray emission is 1.2 x 10~^^ ergs cm~^ s~^ in the 20 — 80 keV band range 
(Fusco-Femiano et al. 2000) and the corresponding luminosity is ~ 10^ ergs s~^. This may 
suggest that the hard X-ray emission observed in clusters so far are not attributed to the 
particles accelerated by turbulence, although our turbulence model is favorable in the view 
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of the maximum energy of electrons (§2.2). In the future, it will be useful to observe the 
spatial distributions of hard X-ray emission in a number of clusters with sensitive detectors. 
If hard X-ray emission is observed from the same region as synchrotron emission but there 
are no shock features such as temperature jumps in this region, it is likely that the emission 
is due to turbulent acceleration of electrons. 

The fluid turbulence could be observed directly by detectors with high spectral reso- 
lution. Figure 3 shows the typical turbulent velocity of the turbulence, Vf, induced by a 
smaller cluster at the distance r from the center of a larger cluster. The maximum velocity 
is ~ 300 km s~^ and it could marginally be detected by observing broadened X-ray emission 
lines with Astro-E2; the turbulence is likely to be observed at the regions where synchrotron 
radio emission is observed. Since it may be difficult to measure the absolute line widths, 
it may be better to compare the line widths in radio-emitting regions which are expected 
to be turbulent with those in other regions of the cluster. Moreover, although the typical 
turbulent velocity is ~ Vt, there should be velocity components between vt and the velocity 
of the smaller cluster, v. These components would be observed as broad wings of the emis- 
sion lines. The spatial scale of the turbulent nonthermal emission region (v^turb) when the 
smaller cluster is at a radius r is shown in Figure 4. 

One concern is whether our assumption of steady turbulence and particle acceleration 
can be justified or not. Figure 5 shows the life span of the turbulence, tturh = ^o/vt- As 
can be seen, ttmb ^ 10^ yr for the outer regions of clusters. On the other hand, we expect 
that the energy distribution of accelerated electrons reaches the steady solution of z/ = 3 in 
the time-scale of ~ = min[ta,te] (§2.1). Since the emission cooling time-scale, tg, does 
not depend on an electron energy spectrum, the relation te < tturb is a sufficient condition 
for the assumption of steady-state. Figure 2 in Sarazin (1999) showed that the electrons 
with 7 ~ 300 have the maximum lifetime of several of 10^ yr, and that electrons with 
7 = p/ (mec) > 10^ have the emission cooling times shorter than 10^ yr. Thus, we think that 
the assumption of steady-state is justified for electrons with 7 > 10^ in the outer region of 
a cluster. Thus, if the radio and hard X-ray nonthermal emission from clusters is due to 
the electrons accelerated by turbulence up to 7 ~ 10^ — 10^, there must be a pool of source 
electrons with 7 ~ 10'^ to be accelerated by turbulence in clusters. These source electrons 
may be provided by the shock acceleration due to cluster mergers and/or gas accretion onto 
clusters that occurred < 10^ yrs ago. In particular, the merger shock wave expected to form 
in front of the merging smaller cluster may provide those electrons. Alternatively, these 
electrons might come from AGNs in clusters, or might be secondary electrons produced by 
the interactions of cosmic-ray ions (Dolag & EnBhn 2000; Blasi 2001; Miniati et al. 2001). 

As the smaller cluster reaches the center of the larger cluster, iturb decreases rapidly. In 
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this region, since the fluid turbulence has not developed fully to accelerate electrons with an 
energy of 7 ~ 10^, the actual nonthermal luminosities of clusters may be smaller than those 
in Figures 1 even if there is a pool of electrons with the energy of 7 ~ 10^. In other words, 
although electrons with higher energy (say ~ 10^) can be accelerated by turbulence because 
of small te, the turbulent acceleration is not a main contributor to the observed nonthermal 
emission. 

Note that all the gas contained in the smaller cluster was stripped before the smaller 
cluster reached the center of the larger cluster in our model. Thus, there is no nonthermal 
emission from the central region of the larger cluster (Figure 1). This and the small time- 
scale of the fluid turbulence in the central region suggest that our model favors peripheral 
radio relics rather than central radio halos. However, the ram-pressure stripping will depend 
on the central structure of the smaller cluster, although the detailed study is beyond the 
scope of this paper. If the smaller cluster has a distinct small scale potential at the center as 
suggested by the observation of Fornax cluster (Ikebe et al. 1996), the gas in it may not be 
removed easily. Also, if the stars of the central galaxy in the small cluster can supply enough 
gas, the gas may survive the ram pressure from the ICM of the larger cluster. Moreover, the 
radiative cooling of the gas in the central region of the smaller cluster may affects the ram- 
pressure stripping (Ritchie & Thomas 2002). In fact, the cluster IE 0657—56 has a rapidly 
moving small component that appears to have passed the cluster center without losing all of 
the gas (Markevitch et al. 2002). 



3.2. Pcirameter Dependence 

Since there are several uncertain parameters in our model, we investigate the dependence 
of the nonthermal emission from merging clusters on these uncertain parameters. 

Figure 6 shows the nonthermal emission when f^, = 0.1 and 0.5. The luminosities are 
small when f^, is small. In particular, the synchrotron luminosity varies more than the inverse 
Compton luminosity. This is because the synchrotron luminosity is affected not only by the 
energy injection from fluid turbulence but also by magnetic flelds. When or vt are small, 
the induced magnetic fields (eq. [60]) is also small, li B < Bcmb, most of the nonthermal 
emission should be the inverse Compton emission rather than synchrotron emission (eqs. [58] 
and [59]). On the other hand, the results do not much depend on //. 

The magnetic fields shown in Figure 2 are roughly consistent with the large values im- 
phed by Faraday rotation measurements (e.g., Lawler & Dennison 1982; Kim et al. 1990; 
Goldshmidt & Rephaeli 1993; Clarke, Kronberg, & Bohringer 2001). However, the magnetic 
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fields must be much smaller than the values implied by the Faraday rotation measurements 
if both nonthermal radio and hard X-ray emission are attributed to the same electron pop- 
ulation with a power-law energy distribution without an upper cutoff (Fusco-Femiano et 
al. 1999, 2000). Thus, we also consider the case where the magnetic fields are small. Fig- 
ure 7 shows the nonthermal emission when b = 0.001. In this case, B is about 30 times 
smaller than that in Figure 2. Although the synchrotron luminosity is not much different 
from that in Figure 1, the inverse Compton luminosity is much larger. This is because the 
total nonthermal luminosity becomes larger (eqs. [13] and [37]) and the fraction of Compton 
luminosity also becomes larger (eq. [59]). Since the maximum inverse Compton luminosity 
is relatively large (~ 10^^-43 gj.g g-i-^^ ^j^g observations of both nonthermal radio and hard 
X-ray emission may be explained by this turbulence resonant acceleration model. If these 
small magnetic fields are correct, then the turbulent region must be different from the region 
where the Faraday rotation is measured. 

In §3.1, we assumed that the fiuid turbulence follows the Kolmogorov law (m = 5/3). 
However, for fully developed MHD turbulence, Iroshnikov (1963) and Kraichnan (1965) 
suggested that the spectral index should be m = 3/2, although this is controversial (Goldreich 

6 Sridhar 1997). On the other hand, if turbulence has not had time to establish a Kolmogorov 
cascade, the index may be m > 5/3. Thus, we discuss the m-dependence of the nonthermal 
emission. In Figure 8, we present the nonthermal emission when m = 1.55 and 1.75. It shows 
that the emission is very sensitive to m. When m — 1.55, the inverse Compton emission 
reaches ~ ig^^"^^ erg s~^, which is consistent with the observed hard X-ray luminosities 
(Fusco-Femiano et al. 1999, 2000) even though the magnetic field is strong {b — 1). (Note 
that if we take m = 1.5, the inverse Compton luminosity is much larger than the observed 
luminosities.) On the other hand, the synchrotron luminosity exceeds the observed values 
of < 10'^^ erg s^^. However, if the maximum energy of electrons is set by fiuid viscosity at 

7 ~ 10*^ (eq. [40]), the observed synchrotron luminosity at ~ 1 GHz may be smaller than that 
shown in Figure 8 while the inverse Compton luminosity is almost the same (Brunetti et al. 
2001a; Petrosian 2001; Fujita & Sarazin 2001). Thus, the model with m < 5/3 can explain 
both the nonthermal radio and inverse Compton emission as a result of particle acceleration 
by turbulence. 

Although the possibility is small (Fujita et al. 2002), a smaller cluster may fall from the 
turnaround radius of a larger cluster (xj = 2). In this case, the smaller cluster has a large 
velocity near the center of the larger cluster compared to the case in §3.1. Figure 9 shows the 
nonthermal emission from the merging clusters. After the smaller cluster enters the virial 
radius of the larger cluster, it produces nonthermal emission. The maximum luminosities are 
not much different from those in Figure 1. Although the smaller cluster has a larger velocity, 
the radius is smaller owing to the larger ram-pressure. However, as shown in Figure 10, the 
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time-scale of the fluid turbulence is small (^turb ^ 10^ yr). This means that the turbulence 
cannot accelerate electrons at 7 ~ 10^ up to ~ 10^, where the observed nonthermal emission 
is produced. Thus, the actual luminosities may even smaller than those shown in Figure 9. 
On the other hand, since the turbulent velocity is large (~ lO'^ km s~^; Figure 11), it would 
be easily observable by high spectral resolution X-ray detectors (e.g., Astro-E2). 



3.3. Preheating 

If the gas distribution of the small cluster has been flattened and the central gas density 
had been decreased by non-gravitational heating before cluster formation (preheating), the 
gas is more easily removed by the ram-pressure when the smaller cluster plunges into the 
larger cluster. Figure 9 shows how the nonthermal emission from merging clusters is affected 
by preheating. Compared with the case of no preheating (Figure 1), the maximum lumi- 
nosities of the nonthermal emission is small especially when the mass of the smaller cluster 
is small. This is because preheating affects less massive clusters significantly and the gas is 
stripped earlier especially when the virial temperature of the smaller of the merging clusters 
approaches 71 (eqs. [49] and [50]). 

Radio observations have shown that the luminosity of diffuse radio emission from clusters 
is a steep function of cluster X-ray temperature or luminosity (Feretti 2000; Liang et al. 2000). 
Since the sample of radio emission is small and less luminous radio clusters may not have 
been observed, this correlation may actiially show that the maximum radio luminosity for a 
given X-ray temperature or luminosity is a steep function of cluster X-ray temperature or 
luminosity. This suggests that the steep correlation may apply only to clusters undergoing 
a major merger (small Ml^^-^,./ Ms^vu)- 

Since the preheating reduces the nonthermal emission from less massive clusters, it may 
help to explain the observed steep correlation. We assume that min(A'/i,vir/^s',vir) ~ 3 and 
it does not depend on M/^ vir- For the no preheating model (Figure 1), the synchrotron 
luminosity from merging clusters of ^ii. = 5 x IQ-'^'^AfQ and Ms^w = 1-5 x IO^'^Mq is 
about 50 times smaller than that from merging clusters of Ml^it = 5 x IQ^^Mq and Ms^u = 
1.5 X IO^^Mq. If the temperature of the larger clusters are represented by Tjcm oc M^^^j^. 
(e.g., Evrard, Metzler, & Navarro 1996), the ratio of the temperature between the two cases 
is 4.6. Thus, the synchrotron luminosity from merging clusters roughly follows Lgyn oc T^^. 
Similarly, for the preheating model (Figure 9) we obtain Lgyn oc Tjcm- Since the virial 
temperature of the smaller cluster of M^^vir = 1-5 x IO^'^Mq is ~ Ti, its merger with the larger 
cluster is significantly affected by the preheating. However, although preheating steepens 
the relation, the power-law index is still smaller than that observed (~ 6; Feretti 2000). 
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3.4. Turbulence Developed in the Whole Cluster 

Numerical simulation done by Roettiger et al. (1999) showed that the fluid turbulence 
in ICM is also developed after a smaller cluster passes the center of a larger cluster. This 
turbulence is pumped by dark matter-driven oscillations in the gravitational potential; the 
merger-induced large scale bulk flows breakdown into turbulent gas motions. In this case, we 
expect the turbulence prevails on a cluster scale, and we can roughly estimate the nonthermal 
emission from the cluster. If we assume that Vt — 1000 km s~^, Iq — 1 Mpc, picM = 
7 X 10-28 g cm-3, TicM = 10^ K, 6 = 1, and m = 5/3, we obtain L^yn = 8 x 10^^ erg s"^ and 
Lie = 1 X 10^^ ergs"^ using the model in §2. Thus, turbulence is strong enough to produce 
the observed synchrotron emission (~ 10'^^ ergs"^), but cannot produce the observed hard 
X-ray emission (~ g-i^ long as m = 5/3 and the magnetic fields are strong 

(6 = 1). However, for the same parameters but with weak magnetic fields (6 = 0.001), we 
obtain Lgyn = 8 x 10^^ erg s~^ and Lie = 9 x 10^^ erg s~^. The luminosities are comparable 
to the observations. Observed central radio halos may be this kind of nonthermal emission 
from the whole cluster. The turbulent velocity, Vt, would be large enough to be detected by 
Astro-E2. 



4. Conclusions 

We have investigated nonthermal emission from electrons accelerated by resonant Alfven 
waves in clusters of galaxies. We assume that the Alfven waves are driven by fluid turbulence 
generated by cluster mergers. We flnd that the resonant Alfven waves can accelerate electrons 
up to 7 ~ 10^; this value is limited by fluid viscosity. Our calculations show that the turbulent 
resonant acceleration can give enough energy to electrons to produce the observed diffuse 
radio emission from clusters if there is a pool of electrons of 7 ~ 10^ in clusters.. On the 
other hand, the observed hard X-ray emission from clusters is explained by the turbulent 
resonant acceleration only when magnetic fields are small (< fiG) or the fiuid turbulence 
spectrum is fiatter than the Kolmogorov law. The fiuid turbulence responsible for the particle 
acceleration would be observed by Astro-E2 in the regions where diffuse radio emission is 
observed. Although non-gravitational heating before cluster formation (preheating) makes 
the relation between radio luminosity and X-ray temperature steeper, our predicted relation 
is still flatter than the observed one. 
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r (Mpc) 

Fig. 1. — Synchrotron (thick hnes) and inverse Compton (thin hnes) luminosities from 
merging chisters when the smaller cluster is at the distance r from the center of the larger 
cluster. The solid, dashed, dot-dashed, and dotted lines correspond to (M^^vir, Ms^w) = 
(5xl0^^ 1.5x10^5), (5xl0^^ 5x10^^), (5xl0^^ 1.5x10^^), (5xl0^^ 5xlO^3)M0, respectively. 
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Fig. 2. — Magnetic fields generated by fluid turbulence when the smaller cluster is at the 
distance r from the center of the larger cluster with the mass of M^^vir = 5 x 10^^ Mq (solid 
hne) and Ml,vit = 5 x 10^^ Mq (dot-dashed hne). 
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Fig. 3. — The velocity of fluid turbulence when the smaller cluster is at the distance r 
from the center of the larger cluster with the mass of Ml^^^ — 5 x 10^^ Mq (solid line) and 
ML,yir = 5 X 10^^ Mq (dot-dashed hne). 
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Fig. 4. — Spatial scale of nonthermal emission when the smaller cluster is at the distance r 
from the center of the larger cluster. Notation is the same as Figure 1. 
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Fig. 5. — Thin lines are the time-scale for decay of fluid turbulence when the smaller cluster 
is at the distance r from the center of the larger cluster. Thick lines are the synchrotron 
luminosities (same as those in Figure 1). Notation is the same as Figure 1. 




Fig. 6.— Same as Figure 1 but for (a) Model Bl (b) Model B2. 
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Fig. 7. — Same as Figure 1 but for Model C. 




Fig. 8.— Same as Figure 1 but for (a) Model Dl (b) Model D2. 
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Fig. 9. — Same as Figure 1 but for Model E {xi = 2). 




Fig. 10.— 



Same as Figure 5 but for Model E (xj = 2). 
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Fig. 11.— Same as Figure 3 but for Model E {xi = 2). 
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Fig. 12. — Same as Figure 1 but for Model F (preheating). 
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Table 1. Model Parameters 



Models 




b 


m 




Preheating 


A 


0.2 


1 


5/3 


0.5 


no 


Bl 


0.1 


1 


5/3 


0.5 


no 


B2 


0.5 


1 


5/3 


0.5 


no 


C 


0.2 


0.001 


5/3 


0.5 


no 


Dl 


0.2 


1 


1.55 


0.5 


no 


D2 


0.2 


1 


1.75 


0.5 


no 


E 


0.2 


1 


5/3 


2 


no 


F 


0.2 


1 


5/3 


0.5 


yes 



